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r  v  r ibution/ 


Abstract 
'  + 

In  this  paper  we  study  a  generalized  second  order  resonance-tunneling 
ordinary  differential  equation,  (z  +  iT)  $"(z)  +  a$'(z)  “  (8l  _  z&22)4>(z)  =  0, 
where  the  parameter  ot  is  real  and  positive,  a  >  0,  and  all  other  parameters 
are  complex.  We  have  obtained  solutions  of  this  equation  when  the  independ¬ 
ent  variable  covers  the  infinite  domain,  -°°  <  z  <  °°,  and  making  use  of  suit¬ 
able  integral  representations  we  have  determined  their  asymptotic  expansions. 
With  the  aid  of  their  leading  terms  the  two  transmission-reflection  problems 
have  been  solved.  Thus,  we  show  that  the  expected  Budden  absorption  occurs 
for  a  =  2N,  N  =  0,  1,  2,....  However,  in  bands  centered  about  a  =  2N  +  1 ,  the 
modulus  of  the  reflection  coefficient  is  generally  larger  than  unity. 

PACS:  03. 40. Kf ,  02.90, +p,  42.10.-s 


1. 


Introduction 


In  a  brief  paper  in  Physical  Review  Letters  1  we  announced  the  preparation 
of  a  lengthier  publication  containing  rigorous  proofs  of  the  outlined  results. 
Thus,  the  present  paper  is  meant  to  serve  the  indicated  purpose.  It  is  well 
known  that  the  propagation  of  waves  through  nonuniform  media  has  been  the 
subject  of  broad  interest  and  study  in  various  areas  of  basic  and  applied 
physics.  We  wish  to  emphasize  the  two  physical  effects  encountered  in  such 
cases;  namely,  the  appearance  of  cutoffs,  i.e.,  points  where  the  local  wave- 
number  vanishes  [k(z)  =  0],  and  resonance  points  where  short  wavelengths 
develop  [(k(z)  <x>] .  In  general,  we  expect  cutoffs  to  result  in  wave  reflec¬ 

tion,  while  resonances  lead  to  wave  absorption.  Because  of  the  seemingly 
opposite  roles  played  by  cutoffs  and  resonances,  it  is  of  considerable 
importance  to  investigate  model  systems  in  which  both  effects  are  present 
simultaneously  and  in  a  sense  compete  with  each  other;  this  is  the  essence  of 
resonance-tunneling  problems.  The  prototype  equation  used  to  approximate 
physical  systems  exhibiting  such  features  is  Budden's  equation.2  This  second 
order  differential  equation  describes  systems  in  which  propagation  regions 
exist  on  both  sides  of  a  resonance,  but  in  which  one  of  the  propagation 
regions  also  contains  a  cutoff.  An  intrinsic  property  of  Budden's  equation  is 
the  fact  that,  for  waves  approaching  the  resonance  (z  =  0)  from  the  cutoff 
side  (z  >  0),  the  modulus  of  the  reflection  coefficient  is  less  than  unity. 
Such  a  result  is  consistent  with  the  intuitive  expectation  than  a  fraction  of 
the  incident  wave  is  absorbed  at  the  resonance,  while  yet  another  portion  can 
tunnel  to  the  propagation  region  on  the  other  side  (z  <  u).  The  converse  sit¬ 
uation  corresponds  to  a  wave  approaching  the  resonance  from  the  left  (z  <  ()) 
and  a  portion  of  which  is  transmitted  to  the  other  side  (z  >  0),  where  the 
cutoff  exists,  while  there  is  no  reflection  at  all.  In  neither  case  is  the 


initial  wave  energy  conserved  in  Budden's  problem  as  well  as  in  our  general¬ 
ized  version. 

One  important  feature  of  Budden's  equation  is  the  absence  of  a  first 
derivative  term.  In  this  paper  we  have  generalized  the  resonance-tunneling 
problem  by  inserting,  in  the  differential  equation,  a  first  derivative 
multiplied  by  an  adjustable  parameter  a  which  is  real  and  positive,  ct  >  0. 

When  we  put  a  =  0  we  recover  Budden's  equation.  We  have  further  generalized 
the  problem  by  allowing  the  presence  of  dissipation,  with  the  result  that  all 
other  parameters  in  the  equation,  as  well  as  the  independent  variable,  become 
complex  numbers.  We  have  solved  the  generalized  equation  in  detail  and 
discovered  that,  as  function  of  the  parameter  ct,  the  reflection  coefficient 
exhibits  an  anomalous  behavior;  namely,  that  the  expected  Budden  absorption 
occurs  for  a  =  2N,  N  =  0,  1,  2,...,  and  that  in  bands  centered  about  a  = 

2N  +  1,  the  modulus  of  the  reflection  coefficient  can  be  larger  than  unity. 

Initially,  our  motivation  for  the  study  of  the  generalized  equation  arose 
from  a  specific  problem3  (electrostatic  whistler  waves  in  magnetized  plasmas 
with  longitudinal  density  gradients),  but  having  encountered  the  extraordinary 
reflection  phenomena  described  above,  we  proceeded  to  carry  out  a  detailed 
analysis  of  the  solutions  of  the  generalized  equation  which  we  present  in  this 
paper.  Thus,  in  Section  2  we  present  the  basic  differential  equation,  and 
then  obtain  a  fundamental  set  of  solutions.  In  Section  3  we  take  up  the  role 
of  Whittaker's  equation  to  study  its  behavior  close  to  the  resonance  (z  =  0) 
and  to  investigate,  tor  a  >  0,  the  small  argument  leading  terms  of  the 
solutions  of  the  generalized  equation.  In  Section  4  we  set  up  integral 
representations  for  the  solutions  of  our  basic  differential  equation.  In 
Section  5  we  compute  and  tabulate  their  asymptotic  leading  terms,  which  we 
need  in  Section  h  to  compute  the  transmission  and  reflection  coefficients. 


Basic  Pi f ferent ial  Equation 


As  explained  in  the  Introduction  we  obtain  our  basic  differential  equa¬ 
tion  by  inserting  in  Budden's  equation  a  first  derivative  term.  Thus,  we  have 

(2  +  iT)  r(z)  +  a<t>'(z)  -  (31  -  2  322)  $(z )  =  U  (2.1) 

which  reproduces  Eq .  (1)  of  Ref.  1.  In  (2.1)  the  parameters  T,  Bj,  and  82 
are  complex  numbers  with  positive  real  parts  and  small  positive  phase  angles, 
which  heur is t ical  ly  account  for  dissipative  processes.  In  the  absense  of 
dissipation,  T  vanishes  and  (Sj,  82)  are  positive  definite  numbers.  The 
parameter  a  in  (2.1)  which  is  the  coefficient  of  the  first  derivative  <j>'(z), 
is  real  and  positive,  0  <  a  <  <»,  and  plays  a  dominant  role  in  the  entire 
development  of  the  present  paper.  Thus,  when  a  =  0  we  recover  Budden's 
equation;  when  a  =  1  we  observe  that  (2.1)  becomes  the  self-adjoint  analogue 
of  the  Budden  problem,  and  when  a  varies  over  the  whole  range,  0  <  a  <  ®,  we 
encounter  the  periodic  variations,  as  functions  of  a,  of  the  reflection  and 
transmission  coefficients  pointed  out  earlier  in  Ref.  1  and  to  be  derived 
rigorously  in  the  sequel.  Mote  that,  as  a  matter  of  notation,  we  employ  in 
(2.1)  the  parameter  T  in  lieu  of  the  parameter  v  appearing  in  Kq .  1  of  Ret.  1. 

The  main  purpose  of  this  paper  is  to  present,  for  the  infinite  domain, 

-®  <  z  <  «>,  the  analytic  and  asymptotic  properties  ot  the  solutions  to  our 
basic  equation.  To  this  end  we  introduce  in  (2.1)  the  following  notation: 

w  =  i£,  %  *  282  (z  +  iT)  ;  (2.2) 


8o  =  S0  +  i  l' 3 2/  2 ,  30  =  81/282  • 


(2.  3) 


We  recall,  from  the  definition  (2.2),  that  in  the  absence  of  dissipation  the 
parameter  82  becomes  positive  definite  and  the  parameter  r  vanishes.  Hence, 
in  this  limit,  5  =  2$2Z  is  a  real  number.  Further  we  recall  that  in  the 

infinite  domain,  <  z  <  “,  we  have  arg  z  =  0  when  z  >  0,  and  arg  z  =  tt  when 

z  <  0,  the  latter  statement  being  a  consequence  of  causality  as  embodied  in 

the  general  definition  £  =  282  (z  +  Thus,  in  the  presence  of  dissipation 

both  82  and  F  are  complex  numbers  with  positive  real  parts  and  (small) 

positive  phase  angles;  that  is  f  =  J  fr| e1  9 1 ,  62  =  |  6^ | e 1  9?- ,  where  0  <  0^,  02  « 

*/2,  from  which  it  follows  that,  in  the  general  case,  arg  £  ♦  02,  (it  +  82) 
as  z  ♦  t®. 

Further  m^re ,  in  what  follows,  we  propose  to  use  w  as  the  independent 
variabLe  and  we  shall  encounter  that  the  multivalued  functions  that  emerge 
will  be  expressed  in  terms  of  ±w  =  ±i£.  Thus,  at  this  juncture  it  becomes  of 
interest  to  ascertain  the  ranges  in  arg(±i£)  which  correspond  to  the  princi¬ 
pal  branch  of  the  functions  in  question.  For  the  purpose  we  assume  initially 
the  absence  of  dissipation,  which  means  that  £  =  282Z  is  real  and  covers  the 
infinite  domain  -<*>  <  £  <  °°,  with  arg  £  =  0  when  £  >  0  and  arg  £  =  it  when  £  <  0 

as  shown  in  Fig.  1,  just  as  in  the  case  of  z  in  the  previous  paragraph. 

First,  to  determine  w  =  i£  we  swing  the  real  axis  in  the  £  plane  through  an 
angle  of  u/2  in  the  counterclockwise  direction  about  the  origin  £  =  0,  which 
makes  w  =  i£  coincide  with  the  axis  of  imaginaries  as  shown  by  the  dashed 
line.  This  procedure  shows  that  the  cut  in  the  complex  £  plane  must  be  drawn 
along  the  negative  imaginary  axis,  and  we  deduce  that  the  principal  branch  for 
a  multivalued  function  of  w  =  i£  is  given  by  the  range 


-n/2  <  arg  (i£)  <  3u/2 


(2.4) 


Then,  to  determine  -w  =  -i £  we  must  now  swing  the  real  axis  in  the  £  plane 


through  an  angle  of  tf/2  in  the  clockwise  direction  about  the  origin  of  coor¬ 
dinates,  which  makes  -w  =  -i £  coincide  with  the  axis  of  imaginaries  as  shown 
by  the  dotted  line.  Hence,  the  cut  in  the  complex  £  plane  must  again  be  drawn 
along  the  negative  axis  of  imaginaries,  exactly  as  before,  except  that  this 
time  we  deduce  that  the  principal  branch  of  a  multivalued  function  of  -w  =  -i £ 
is  given  by  the  range 

-ir/2  <  arg  (-i£)  <  3ir/2  ,  (2.5) 

which  should  be  contrasted  with  (2.4).  We  observe  in  Fig,  1  that  the  cut 
along  the  negative  axis  of  imaginaries  separates  the  dashed  and  dotted 
segments  in  the  lower  half  plane  of  w  =  i£  and  -w  =  ~i£,  respectively,  and  we 
note  that  the  corresponding  segments  in  the  upper  half  plane  actually  coin¬ 
cide  with  each  other. 

In  the  presence  of  dissipation,  as  we  have  already  seen,  we  have 
£  =  202(z  +  iO,  where  &2  and  r  are  complex  parameters,  and  asymptotically 
(|z|  °°)  this  definition  becomes  £  ~  202z.  Confining  our  attention  to 

positive  values  of  z,  0  <  z  <  °°,  we  compute  the  complex  numbers  w  =  ±i£,  which 
are  shown  in  Fig.  I  by  two  black  dots  equidistant  from  the  origin  and  lying 
on  a  straight  line  inclined  at  an  angle  02  =  arg  S>2  with  respect  to  the  axis 
of  imaginaries.  We  observe  that  the  transition  from  w  in  the  upper  half  plane 
to  -w  in  the  lower  half  plane,  or  vice  versa,  is  achieved  by  the  semicircle 
drawn  to  the  right,  which  conforms  to  the  convention 

(-v)a  =  wa  elws.a)  e  =  sgn  Re{w},  (2.5) 

where  a  is  an  arbitrary  coefficient.  In  Fig.  ,  the  point  w  is  such  that 
Re{w}  <  0  and  Im{w}  >  U;  that  is,  w  is  a  positive  pure  imaginary  (02  =  IJ)  or 


else  lies  on  the  second  quadrant  of  the  complex  ^  plane,  with  corresponding 
statements  for  the  point  -w  on  the  lower  half  plane. 

Thus,  in  making  use  of  w  =  i £  as  the  independent  variable  in  (2.1),  we 
obtain,  writing  Y(w)  instead  of  t>(z),  the  generalized  equation 

w  Y" (w )  +  aY'(w)  +  (i30  -  w/4 )  Y(w)  =  0  ,  (2.7) 

which  we  proceed  to  solve  in  terms  of  confluent  hypergeometric  functions. 
Putting  Y(w)  =  e-w/2y(w),  it  is  found  from  (2.7)  that  y(w)  satisfies  the 
canonical  form  of  Hummer's  equation 

w  y" (w)  +  (b  -  w)  y'(w)  -  ay(w)  =  0  (2.8) 

with  parameters 

a  =  a/2  -  iB0,  b  =  a,  b  -  a  =  a/2  +  i(30  .  (2.9) 

Two  linearly  independent  solutions  of  (2.8),  valid  for  all  possible  values  of 
the  parameters  (a,b)  are4 

y5(w)  =  U( a ,  b;  w) ,  y?(w)  =  e«  U(b  -  a,  b;  -w),  (2.10) 

in  which  U(a,  b;  w)  is  a  multivalued  solution  to  Hummer's  equation  for  which 
the  principal  branch  is  usually  taken  as  -it  <  arg  w  <  it.  Consequently,  a 
fundamental  set  of  solutions  of  (2.7)  is  given  by 


Y_(  w )  =  -w  /  2  y  -  (  w  )  =  l--w  2  l'(  a  ,  b  ;  w)  =  e  w/2  U  (  at/  2  -  i  Bc  ,  a;  w)  ;  (2.11) 

Y+(w)  =  e-w/-y7(w'l  =  ew  2  L'(  o  -  a,b;  -  w)  =  ew/2  l' ( a/ 2  +  i80,  a;  -  w),  (2.12) 

Hi  which  the  latter  forms  apply  specifically  to  the  parameters  defined  by 
(2.9).  It  is  noteworthy  to  point  out  that  these  parameters  (2.9),  and  the 
corresponding  solutions  (2.11)  and  (2.12)  satisfy  the  transformation 

a  «•-*•  b  -  a,  b  b,  w  *■  +  -w;  (2.13) 

that  is,  Y_(w)  as  given  by  (2.11)  can  be  turned  into  Y+(w)  as  given  by  (2.12), 

or  vice  versa,  by  merely  applying  the  transformation  (2.13). 

For  brevity  in  writing,  it  is  convenient  to  introduce  from  (2.9)  the 
parameters 

a±  =  a/2  -  i3o,  b  =  a,  (2.14) 

with  the  aid  of  which  the  solutions  (2.11)  and  (2.12)  become 


Y-(v) 

=  e-w/2  LI ( a _ ,  a;  w) 

|arg  (w)j  <  it; 

(2. 

13) 

Y+(w) 

=  ew/2  L!(a  +  ,  a;  -w) 

I  arg  (-w)j  <  it. 

(2. 

16) 

which  the  reader  will  recognize  in  Kef.  1  as  Kqs.  (5)  and  (4),  respectively. 

It  is  readily  verified  that  the  parameters  (2.14)  and  the  solutions  (2.15)  and 
(2.1b)  satisfy  the  transformation  (2.13),  i.e.,  a+  a_,  a  a,  w  <-♦  -w , 
which,  with  w  =  i%,  implies  that  (2.13)  can  be  condensed  into  a  single 
reversible  transformation  i  ►  -l  ,  wherever  i  appears  explicitly. 


—  H  - 

3 .  The  Role  ot  Whittaker's  liquation 

In  Che  sequel  we  shall  often  find  it  more  convenient  to  deal  with  the 
complex  independent  vat  table  2  =  26n(z  +  ll),  as  defined  by  (2.2)  than  with 
the  independent  variable  w  =  i; ,  Accordingly,  writing  Z(^)  instead  of  Y ( 1  % ) 
in  (2.7)  yields 

+  aZ  ’  (  r, )  +  (5/4  -  30)  Z(  5)  =0.  (3.1) 

We  can  show  that  (3.1)  is  related  to  Whittaker's  equation  by  putting 

Z(  U  =  cr!X'1  u(P),  (3.2) 

where  u(5)  satisfies  the  Helmholtz  equation 

u"(y)  +  Q(5)  u (5)  =  0,  (3.3) 

in  which 

<K5)  =  1M  "  Sc/r,  +  a(2  -  *)Mr.2  (3.4) 

plays  the  role  of  an  effective  uotential  and  is  real  in  the  absence  ot  dissi¬ 
pation;  that  is,  when  T  =  0  in  (2.3).  It  is  clear  from  (3.4)  that,  when  a  =  U 
or  a  =  2,  the  two  cases  exhibit  identical  potentials  dominated  in  the  vicinity 
ot  the  resonance  (5  =  0)  by  the  anti-symmetric  term  c~ 1 .  On  the  other  hand, 
when  a  differs  from  these  special  values,  a  significant  change  develops  near 
the  resonance.  This  time,  for  0  <  a  <  2,  it  acquires  a  symmetric  term  that 
behave.-,  as  and  for  2  <  a  <  00  the  symmetric  term  becomes  .  It  appears 

that  the  presence  ot  the  tirst  derivative  in  (3.1)  is  responsible  tor  this 
modification  ot  the  potential  0(?)  In  the  vicinity  of  the  resonance  and  is  the 
principal  reason  tor  the  anomalous  reflection  properties  mentioned  in  Kef.  1 
arid  m  the  Introduction,  prouerties  that  we  demonstrate  rigorously  in  the 


s  enue 1 . 


-9  - 

■  juite  generally  we  observe  chat,  when  ct  =  0,2,  Che  equation  t}(£)  =  U  has 
only  one  root,  namely  ^  +  =  -»80,  which  is  the  classical  cutoft  in  Budden's 
equation.  When  a  is  different  from  these  two  special  values,  the  equation 
o( )  =  '  is  a  quadratic  and  has  two  distinct  roots,  and  f+,  the  intercepts 

of  the  corresponding  curves.  It  we  assume  that  j  80|  is  large  or,  more 
precisely,  that  | A  3  3  2  j  >>  ja(  a  -  2);,  we  deduce  that  the  intercepts  are  given 
approximately  by 

■>-  3  _a( a  -2)/4S0;  3  98o  +  a(  a  -  2)/4So>  (3.5) 

which  satisfy  the  general  conditions  for  the  roots  of  a  quadratic, 

’-+?+=  ^B0,  £_  *  U  =  a(  2  -  a),  (3.6) 

and  which  teLl  us  that,  under  the  present  assumptions,  the  cutoff  point  at 
*>+  =  430  (for  a  =  0,2)  is  shitted  only  by  a  relatively  small  amount  as  a  is 
changed.  Conversely,  we  notice  from  (3.5)  that  the  cutoff  is  relatively 
small  but  greatly  dependent  on  ou  Finally,  we  observe  that,  when  a  =  1,  which 
corresponds  to  the  self-adjoint  analogue  of  the  Budden  problem,  acquires 
its  maximum  value,  l  /  4  S0 ,  while  z,+  its  mininum,  4  80  -  1/4  S0 ,  statements  which 
imply  that  3, ,  is  assumed  real. 

To  illustrate  the  essential  features  discussed  in  the  preceding  para¬ 
graphs  we  present  in  Fig.  2  the  curves  for  0(F),  as  given  by  (3.4),  assuming 
for  the  moment  that  S0  is  real.  In  this  way  we  obtain,  putting  80  3  1  for 

convenience,  the  curves  corresponding  to  a  =  0,2,  a  =  1,  and  a  =  3.  We  observe 
that  ttie  transition  from  the  case  a  =  0,2,  to  either  the  intermediate  region 
;)  <  a  <  2  or  to  the  extended  domain  l  <•.  a  <  ”,  is  non-uniform.  Final  Ly,  we 
wish  to  recall  that  the  mam  thesis  of  the  present  paper  is  the  enhanced 


reflection  that  we  observe  in  certain  periodic  bands  in  a,  in  which  it  appears 


that  the  vicinity  of  the  resonance  at  5  =  0  is  stimulated  to  emit  waves 
spontaneously.  To  the  best  of  our  knowledge  there  is  no  intuitive  way  of 
extracting  this  periodic  behavior  from  a  study  of  Q(  5)  from  the  curves  of 
Fig.  2  or  analytically  from  the  definition  (3. A). 

Having  studied  the  behavior  of  the  function  Q(5)>  it  becomes  of  interest 
to  examine  the  small  argument  leading  terms  of  the  solutions  (2.13)  and  (2.16) 
which  we  rewrite  as 

Y±(i£)  =  e-i  ’/2U(a  +  ,  a;  ±U),  (3.7) 

where,  it  is  recalled,  a±  =  a/2  ±  i30  and  i£  =  2i62(z  +  i?)  =  2i  S2Z  -  202?, 
with  j  2 S2 T |  <<  1.  The  exponentials  in  (3.8)  thus  become 


±15/2 
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iiS^z  -t-BoT  T 

e  -  •  e  1  - >  e  1 

z  0 


1+82 r  +  ..., 


(3.8) 


which,  because  |  02^’ I  «  3,  may  be  equated  to  unity  in  a  small  argument  (z  ♦  U) 
expansion.  Consequently,  we  need  only  look  at  the  Kummer  functions  in  (3.7). 
Thus,  when  a  is  not  zero  or  an  integer  (i.e.,  a  *  0,1,2,...),  we  express 
U(a±,  a;  ±i5)  as  a  linear  combination  of  Kummer  functions,6  namely: 


U(a+,  a;  ±i5)  = 


rq  -  a) 

T(a±  -  a  +  1) 


M(af,  a;  +i  5) 


+  —ry  "c1--  (^15)1~a  M(a±-a-rl ,  2-a;  ±15).  (3.9) 

1  (a+ ) 

To  obtain  the  small  agrument  expansions  for  these  values  of  a  we  replace  the 
Kummer  functions  on  the  right  of  (3.9)  by  unity,  which  is  their  leading  term, 


V  V 


to  obtain 


l'(a±,  a;  *i '  ) 


•  v  i  ZJzI  +  1  v  a  ~  1  > 
r(a+  —  oi+l)  T(a+) 


(3.  10) 


from  which  we  deduce, 


for  0  <  a  <  1 : 


Y±(i^) 


„  £U  -  a) 

r(  a  +  -a  +  1)  ’ 


( 3.  1  l) 


for  1  <  a  < 


Y±(ir,) 


"(3  -  O 

T(a±) 


(HUl"a  , 


(3.  Id) 


which  show  that  the  leading  term  (3.11)  is  a  constant  (independent  of  ■" ) , 
whereas  the  terms  (3.12)  behave  as  reciprocal  powers  of  (i^). 

When  a  =  0,1,2,...,  the  formula  (3.9)  diverges,  and  we  must  resort  to  a 
limiting  process  that  leads  to  the  so-called  logarithmic  case.*1  Thus,  it  we 
put  a  =  n  +  1,  n  =  1,2,...,  whence  a*  =  (n  +  l)  ±  i 30 ,  the  small  argument 
leading  term  becomes 

U(a  +  ,  n  +  1;  ?  ir,)  =■  -^r - ^  (  n^)-0  ,  (  3.  13) 

1  (ae) 

from  which  we  deduce,  for 

a  =  2,  n  =  1,  a±  =  1  ±  i80;  Y±(icJ  *  r(ae)~  »  (3.U) 

in  full  agreement  with  (3.12).  However,  the  same  formula  yields 

a  =  1,  n  =  0,  a+  =  t  iSQ;  Yt(i^)  a  -  j^TT)  ln  ( ?i  •»)  (3.  15) 

a  =  0,  n  =  -1,  a+  =  ±  i 30 ;  Y±(i')  51  T^aTJ  - 1  +  ^ ^  •  (3.  Id) 

which  are  the  logarithmic  results  mentioned  earlier  that  cannot  be  obtained 
from  (3.12)  and  do  not  blend  smoothly  with  the  values  obtained  tor  the  range 


0  <  a  <  1 


To  obtain  integral  representations-'  of  the  solutions,  we  introduce  the 
Laplace  kernel  K(?,t)  =  e*’11  and  apply  to  (3.1)  the  standard  procedures'’  that 
abide  when  all  the  coefficients  of  the  various  terms  ot  the  differential 
equation  are  at  most  linear  tunctions  of  the  independent  variable.  Thus,  we 
assert  that  the  integral  representations  of  the  solutions  ot  (3.1)  may  be 
written  as 


Z(£)  = 


(t  -  |) 


(4.1) 


where  A  is  an  arbitrary  normalization  constant  to  be  determined  a  posteriori, 
a+  and  a_  are  the  parameters  of  the  problem  as  defined  by  (2.14),  and  C  is  a 
suitable  path  of  integration  in  the  complex  t  plane,  which  must  be  chosen  to 
conform  with  the  bilinear  concomitant 


PU.O 


(t  +  y) 


a_  i  £t 
e 


(4.2) 


That  is,  the  contour  C  must  be  such  that  P(S,t)  returns  to  its  primitive  value 
after  t  has  described  a  closed  circuit,  or  else  the  path  of  integration  is 
drawn  between  definite  limits  such  that  P(S,t)  vanishes  at  each  limit.  In  the 
sequel,  we  assume  that  the  cuts  in  the  complex  t  plane  are  chosen  along  the 
real  axis;  they  extend  from  t  =  1/2  to  00 ,  and  from  t  =  -1/2  to  -°°.  The 
branch  points  at  t  =  —1/2  cannot  be  connected  to  each  other  by  a  single  branch 
cut  because  the  coefficients  of  the  two  binomials  in  the  integral  (4.1)  are 
different  from  each  other.  The  reader  wi L 1  readily  recognize  that  (4.1)  and 
(4.2)  reproduce  Eqs.  (8)  and  (9)  of  Ref.  1. 

We  have  adopted  for  convenience  the  latter  procedure  (the  use  ot  fixed 
limits  of  integration)  noting  that  the  bilinear  concomitant  (4.2)  vanishes  at 


I* 


4 


3 


the  branch  points  t  =  ±1/2,  provided  we  have  Re { a  + }  >  0,  which  excludes 
Budden's  problem  (a  =  0)  that  we  treat  properly  later.  We  also  note  that 
P  +  0  as  | t |  *  “  along  a  ray  inclined  (with  respect  to  the  axis  of  reals)  at 
an  angle  (y  -  <£),  where  j>  satisfies  the  condition8 


■j  +  <  arg  £  <  y  + 


(4*3) 


that  guarantees  the  convergence  of  the  integral  (4.1)  as  |t|  grows  without 

limit;  for  example,  when  $=  0,  we  have  -  y  <  arg  £  <  y,  which  includes 

IT  3  IT 

arg  £  =  0,  and  when  <p  =  n,  we  obtain  y  <  arg  £  <  — y,  which  includes  arg  £  =  it. 
In  fact,  in  the  bilinear  concomitant  (4.2),  we  choose  arg  {i £t }  =  tt ,  which  is 
tantamount  to  choosing  the  path  of  steepest  descents,  from  which  we  deduce 
asymptotically,  with  arg  t  ~  y  -  that  =  arg  £.  Likewise,  we  observe  that 

(4.3),  provided  only  that  the  range  of  <j>  is  properly  chosen  as  shown  below, 
guarantees  the  basic  condition 


0  £  arg  £  £  tt  , 


(4.4) 


which  follows  from  (2.2)  and  subsequent  discussion  in  the  limit  of  vanishing 
dissipation,  i.e.,  £  =  2B2Z  is  real. 

Accordingly,  we  define  the  function  Z+( £)  by  means  of  the  rectilinear 
integral  representation 


Z.(£)  =  A  / 

1/2 


l/2+i“e 


i£t  1  a  -1  1  a 

e1^  (t  -  y)  +  (t  +  y)  “  dt  , 


(4.5) 


which  starts  at  t  =  y  and  proceeds  to  infinity  along  a  ray  inclined  at  an 
angle  with  respect  to  the  axis  of  imaginaries,  where  satisfies  the  range 


-14- 


7T 

~  J  <  <P  < 


3n 

2  * 


(4.6) 


which  contains  $  =  arg  £  =  0  and  $  =  arg  £  =  it  in  accordance  with  the  require¬ 
ment  (4.4).  In  addition,  we  observe  from  (4.6)  that,  in  varying  $  about  t  =  y, 

the  rectilinear  path  of  integration  does  not  encroach  upon  the  other  branch 
point  at  t  =  -  y.  If  <(>  were  to  be  increased  beyond  the  limits  stated,  it  would 
become  necessary  to  make  the  path  of  integration  circle  the  branch  point  at 
t  =  -  y,  with  an  additional  contribution  to  the  answer  which  does  not  belong 
in  (4.5).  Then,  if  we  substitute  the  specified  limits  (4.6)  into  the  values 
of  <J>  in  (4.3)  we  find  for  arg  £  the  permissible  range  -u  <  arg  £  <  2ir,  which 
contains  symmetrically  the  basic  condition  (4.4).  Finally,  subtracting  y 
from  each  term  of  the  preceding  set  of  inequalities  we  obtain  the  result 
-  y~  <  ar8  (-i?)  <  y— ,  which  comprises  symmetrically  the  customary  principal 
branch  -n  <  arg  (-i£)  <_  tt  for  the  function  U(a+,  a;  -i£)  that  we  associate 
with  (2.16)  and  (4.5).  However,  as  demonstrated  in  Fig.  1  and  accompanying 
discussion,  the  argument  range  for  the  principal  branch  of  multivalued 
functions  of  (-i£)  is  given  by  (2.5);  that  is,  -n/2  £  arg  (-i£)  <  3 tt/ 2 ,  which 

is  the  range  we  shall  adopt  henceforth  for  the  Kummer  function  U(a  +  ,  a;  -i£), 

and  which  contains  the  range  -it/ 2  arg  (-i£)  n/2  that  corresponds  to 

(4.4). 


To  evaluate  the  integral  (4.5),  we  change  the  variable  of  integration 
from  t  to  u  by  putting 

i£t  =  i £ ( t  -  y  +  y)  =  i £/ 2  +  i£  (t  -  =  i£/2  -  u, 


where  -u  =  ei^u  ; 


(4.7) 


elTtu  =  i  £(  t  -  y),  t  -  y  =  e^u/K;  t  +  y 


u 

i  £  ’ 


1 


dt  =  (ei1T/i  C )du ; 


(4.8) 
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0;  t 
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i°°e 


~iP 
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i  8 
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»Uk“arg  r-  ‘  *! 


(4.9) 


which  states  that  the  upper  Limit  in  (4.5)  tends  to  0Qei-6>  where  B  =  arg£  -  <j> 
and,  according  to  (4.3),  3  is  an  acute  angle,  positive  or  negative. 
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(  4 .  1  U  ) 


To  proceed  with  the  evaluation  of  the  integral,  we  insert  (4.8)  and  (4.9)  into 
(4.5)  to  obtain 


z+(5) 


A  (eiff/i^)a+  e1 


;/2 


°°e 


10 


-u  a+-' 
e  u  T 


(1  ~ 


u  a."1 


ft> 


du 


=  Af(a+ 


)e 


iC/2 


U(a  +  ,a;-iC) 


~  Ar(a+)elC/2(-i?)  a+  2F0(a+,  1  -  a_;  1^) ,  (4.11) 

where  the  result  in  the  second  line  stems  directly  from  (2.16)  and  states  that 
Z+(C)  in  (4.11)  is  proportional  to  Y+(lO;  that  is,  Z+(  £ )  =  AT(a  +  )  Y+(iO, 
where  the  constant  of  proportionality,  AT(a+),  arises  from  the  integral 
representation  itself.  The  last  form  in  (4.11)  is  deduced  from  the  first 
line,  where  the  integral  representation  makes  use  of  u  in  (4.7)  as  the 
variable  of  integration.  Thus,  expanding  the  binominal  (1  -  u/iO3""*  into  a 
partial  series  plus  a  remainder;  interchanging  the  order  of  integration  and 
summation,  and  integrating  term  by  term,  we  obtain,  according  to  Watson's 
Lemma, ’  an  asymptotic  series  in  the  sense  of  Poincare,  which  is  condensed  in 
the  last  form  of  (4.11),  Here,  the  function  2F()  is  a  generalized  hyper¬ 
geometric  series  written  in  Pochhamer's  notation  as  modified  by  Barnes.10  It 


/.v'!-/.-,'0'.- ..... 


A  mK  ktak  ^ 


vL‘4l 


represents  symbolically  an  asymptotic  series  in  reciprocal  powers  of  (l?)  that 


must  be  terminated  after  a  finite  number  ot  terms  tol Lowed  by  a  remainder.  It 
is  recalled  that  in  such  series,  tile  remainder  is  ot  the  order  ot  the  first 
discarded  term,  and  thaL  the  functions  a*-, )  have  unity  as  tile  Leading  term. 
Comparing  the  last  two  forms  in  (4.11)  Leads  to  the  asymptotic  relation 


L'  ( a 


«;  -  irJ 


(-lrJ 


(-4.  1  2) 


which  apparent ty  was  first  noted  by  fcirdelvi  in  his  definition  of  the 
corresponding  Kummer  function.11 


Similarly,  the  companion  function  2_(-)  is  defined  by  means  of  the 

! 

rectilinear  integral  representation 


2  (4) 


-1/2-t-i-e  1  f 

f  e1,C 

-1/2 


4>a*‘ 


(t 


a  -1 

) 


dt 


(4.13) 


which  starts  at  t  =  -1/2  and  proceeds  to  infinity,  as  before,  along  a  rav 
inclined  at  an  angle  -f>  with  respect  to  the  axis  ot  imagiuaries,  where  f> 
satisfies  this  time  the  range 


-it  <  i>  <_  ^  ,  (4.14) 

which  has  been  so  chosen  because  it  contains  f>  =  arg  '  =  U  and  ti  =  arg  r,  =  n 
in  accordance  with  the  requirement  (4.4).  However,  this  time,  when  the 
rectilinear  path  ot  integration  is  rotated  about  t  =  -1/2,  say  f  rum  t>  =■  t  to 
$  =  0,  it  becomes  necessary  to  make  the  path  ot  integration  circle  the  branch 
point  at  t  =  1/2,  giving  rise  to  an  additional  contribution  which  is  not 

contained  in  (4.1))  and  which  is  part  ot  the  process  >1  analytic  continuation, 
a  matter  to  which  we  return  in  Section  a  when  wr  uldiess  the  t r anstn ss  i  m- 
reflection  problems  that  we  seek  to  solve.  Proceeding  as  belole,  II  We 


substitute  the  specified  limits  (4.14)  into  the  values  of  in  (4.3),  we  find 


for  arg  Z,  the  permissible  range  — 3n/ 2  <  arg  Z  £  3m/2,  which  contains  the 
basic  condition  (4.4).  Finally,  adding  n/2  to  each  term  of  the  preceeding  set 
of  inequalities,  we  obtain  -t  <  arg  (id  <  2n,  which  contains  the  principal 
branch  -  y  <  arg  (id  £  for  the  function  l'(a_,  i;  id  that  we  associate  with 

T 

(2.12)  and  (4.13).  Notice  that  this  range  contains  —  arg  (ir,)  ,  which 

corresponds  to  (4.4),  and  that  it  agrees  with  (2. a),  as  demonstrated  in  Fig.  1 
and  accompanying  discussion. 


To  evaluate  the  integral  (4.13)  as  before,  we  change  the  variable  of 
integration  from  t  to  u,  but  this  time  we  put 


i  %t  =  ir,(t  +  y  -  y)  =  -i  f/2  +  iZ(i  +7)  =  -i  CJ1  -  u, 
i  " 

where  -u  =  e  u; 


(4.15) 


i  it  , ,  1  1  in, 

e  u  =  i-,(t  +  y ) ,  t  +  ~2  ~  e  u / i  ;  t 


d 


dt  =  (ein/irjdu; 


(4.  16) 


t  =  -  7,  u  =  0;  t  = 


1  ,-it>  ,  , 

-  7  +  i°°e  ,  u  =  i  u  |  e 
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ri  i  (arg  Z~t>) 
;  e 


(4.17) 


which  implies  that  the  upper  limit  in  (4.13)  tends  to  asef°,  where  B  =  arg  Z  ~ 
is  an  acute  angle,  positive  or  negative.  Thus,  proceeding  as  before,  we  have 
f  rom  ( 4 . 13) 
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=  Aein(a_1) 
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.  i  t  (  a- 1 )  «  -iZl'2 , '  „  -a-  „ 
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where  the  expression  in  the  second  line  stems  directly  from  (2.15)  and  states 


that  Z _(£)  in  (4.18)  is  proportional  to  Y_(i£),  that  is,  Z_( F) 


Ae^7^01  1 )  T(a_)  Y_(  i  5  )  ,  where  the  constant  ot  proportionality,  Ae1  a~l)  T(a_)  , 


arises  from  the  integral  representation  itself.  The  last  form  in  (4.18)  is 


deduced  from  the  first  line  employing  the  same  methods  that  led  to  (4.11), 


Here,  again,  the  function  2f(j  ls  an  asymptotic  series  in  reciprocal  powers  of 


(i?)»  whose  leading  term  is  unity.  Finally,  comparing  the  last  two  results  in 


(4.18)  leads  to  the  asymptotic  relation 


U(a_,  a;  H)  -  (H)  3“  2FQ(a_,  1  -  a  +  ,  -  jj)  , 


(4.19) 


which  should  he  contrasted  with  (4.12),  and  which  again  is  due  to  Erdelyi's 


formulation. 


To  determine  the  constant  A  appearing  in  (4.11)  and  (4.18),  we  demand 


that  the  functions  Z+(£)  and  Z_(£)  satisfy  the  reversible  transformation 


(2.13).  That  is,  in  accordance  with  the  last  sentence  of  Section  2,  we  demand 


that  Z+(f)  =  Z_(F)  and  Z_(5)  =  Z+(£),  where  a  bar  over  a  symbol  signifies  that 


we  have  imposed  the  reversible  transformation  i  *■*  -i ,  wherever  i  appears 


explicitly.  Thus,  since  the  normalization  constant  A  is  arbitrary,  we  put 


A  =  e 


-in(a-l)/2 


(4.20) 


which,  when  inserted  into  (4.11)  and  (4.18),  yields  the  results 


-i w( a  -  1 ) / 2  .  i?/2 

/.  +  (0  =  e  F(a  +  )e  h(a  +  ,  a;  -i£) 


— i  TT  ( ot—  1 )  /  2  „  i^/2  -a+  1 

e  r(a+)e  (-K)  2F0(a  +  ,  1  -  a_;  _)  ; 


(4.21) 
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e  r(a_)e  U?)  2F()(a_,  1  ”  a+;  “IT  ' 


(4.22) 


chat  clearly  satisfy  the  reversible  transformation  i  «-♦  -i . 

When  a  =  0  in  (2.L),  we  recover  Budden’s  equation,  and  the  parameters  of 
the  problem  become  a+  =  ±  ig  ,  which  means  that  the  integral  representations 
(4.5)  and  (4.13)  cannot  be  used  because  the  restriction  Re{a  +  }  >  0  is  being 
violated.  To  obviate  this  difficulty,  we  make  use  of  the  identity^ 


U(a,0;w)  =  wU( 1  +  a,2;w), 


(4.23) 


which  allows  us  to  write  down  the  required  solutions  for  a  =  0  with  the  aid  of 
(2.16)  and  (2.15),  respectively.  Thus,  we  have 


Z+(C)  =  -  i T ( 1  +  i8o)e1?/2U(iBo,  0;  -i £) 


=  — i T (  1  +  ii  )el5/2(-i?)U(  1  +  il,2;-i?) 


-i  T(  1  +  iS  )elC/2(-i^)  l8o  ,Frt(l  +  iB_,  if?.;  -pr)  ; 


(4.24) 
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Z  (?)  =  iT(  1  -  ii  )e_lC/2U(-ii  ,  0;  U) 


i r ( 1  -  iSo)e_U/2(iC)U(l  -  iB0,  2;  if) 


ir(i  -  ii0)e-i5/2(iC)i8o  2F0(1  -  u0,  -i 80 ;  -jj)  ; 


(4.25) 


k' 


which  should  be  contrasted  with  their  counterparts  (4.21)  and  (4.22)  valid  for 

* 

a  >  0,  and  in  which  we  have  made  use  of  the  Erdelyi  asymptotic  relations 
(4.12)  and  (4.19),  respectively. 


/  ■ 


The  most  important  results  of  the  preceeding  paragraphs  are  illustrated 
graphically  in  Fig.  3.  Thus,  Figs.  3(a)  and  3(b)  depict  the  t  plane  and 
emphasize  the  rectilinear  half-ray  paths  of  integration  for  the  functions 
Z+(f)  an  Z_(^),  respectively,  which  we  have  chosen  to  draw,  in  the  absence  of 
dissipation,  parallel  to  the  axis  of  imaginaries.  That  is,  making  sure  that 
in  both  cases  the  causality  condition  (4.4)  is  preserved.  This  is 
accomplished  by  exhibiting  the  points  f  =  arg  c,  =  0  and  <p  =  arg  ?  =  if  on  the 
paths  of  integration,  and  joining  them  by  semicircles  drawn  to  the  right  of 
the  vertical  axis  as  part  of  the  unit  circles  which  correspond,  respectively, 
to  the  limits  (4.b)  and  (a. 14).  Figs.  3(c)  and  3(d)  represent, 
respectively,  the  full  argument  ranges  for  (-i*j)  and  (i%),  which  are  the 
arguments  appearing  in  Z  +  (  f, )  and  Z_(£).  In  addition,  the  portions  of  the 
spirals  drawn  with  t  tiles  lines  represent  the  respective  ranges  corresponding 
t  )  (4. a),  i.e.,  1  <  arg  r,  \  t  .  In  the  presence  of  dissipation,  as  shown  m 


i 


(2.2)  and  subsequent  discussion,  the  basic  argument  range  (4.4)  becomes,  as 
z  grows  without  limit,  3a  <•  r  <  +  it,  where  On  =  arg  32*  The  diagrams 


5 .  Asymptotic  Leading  Terms 


In  this  section  we  are  concerned  with  the  tabulation  of  the  Leading  term 
asymptotic  expansions,  expressed  as  functions  of  i £  =  2i82(z  +  12)  and 
eventually  as  functions  of  the  physical  independent  variable  z,  -®  <  z  <  «>,  a 
was  done  in  Ref.  1.  The  former  objective  is  achieved  by  replacing  the 
functions  2^0  by  unity.  Thus,  from  (4.21)  and  (4.22),  we  have  the  asymptotic 
leading  terms  with  simplified  normalization  factors: 


Z  +  (rJ  -  i(  1"j)r(a  +  )e1’/^(-i()"a+  ; 


(r)  ^-i(1-a)r(a  )e'lU2(il)'a“  , 


(5.1) 

(5.2) 


which  are  valid  for  the  parameters  a  *  =  -j  ±  i 3q  and  a  >  U.  In  the  case  of 
i  =  Hadden's  problem,  we  obtain  the  asymptotic  leading  terms 

Z  (f.)  -  -12(1  +  IS  )el5/2(-iO“i8o  ;  (3.3) 

+  o 

7.  (  rJ  '12(1  -  iS  )e_U/2(H)iS°  .  (5.4) 

o 

In  both  pairs  of  equations  the  adopted  notation  means  that  the  subscript  plus 
in  Z+(  i)  is  associated  with  progressive  waves,  whereas  the  subscript  mi  mis  in 
t.  -  (  f, )  with  regressive  waves.  This  wave  character  of  the  preceeding  solutions 
is  evident  asymptotically;  i.e.,  for  |i5|  >>  1,  where,  it  is  recalled, 

\  =  2 3 2 ( z  +  i r )  in  which  82  and  T  are  complex  parameters.  In  the  absence  ot 
dissipation,  82  is  real  and  positive  and  T  vanishes;  therefore  in  this  case. 


A  is  real. 


First,  we  examine  to  what  extent  the  presence  or  dissipation  afreets  the 
asymptotic  leading  terms  listed  before,  especially  as  regards  the  powers  of 
(i£).  For  the  purpose,  we  assume  temporarily  that  z  =  jzj ,  which  is  positive 
definite,  and  define  from  (2,2) 


which  makes  s  a  complex  number  with  positive  real  part  and  small  positive 
phase  angle;  i.e.,  Im{iSi  >  U  because  arg  %  =  arg  82.  With  this  understand¬ 
ing,  e.g.,  the  computation  of  the  power  (ir,)~a,  where  a  is  an  arbitrary 
exponent,  real  or  complex,  becomes 

(U)_a  =  (H  -  n)“a  =  (K)~a  (L  -  ^)_a 

i? 

=  (i?)"a  rl  +  —  -  (U)"a  (5.6) 

H 

which  is  valid  to  lowest  order  in  reciprocal  powers  of  (is)  because  n  =  2B^V 
is  a  constant.  Finally,  according  to  the  definitions  (2.2)  and  (5.5),  when 
z  >  U,  we  have  ^  =  s,  aind  when  z  <  0,  then  5  =  -£• 

Next,  we  recall  that  the  confluent  hypergeometric  functions  employed  in 
(4.21)  and  (4.22)  for  a  >  0,  as  well  as  their  counterparts  (4.24)  and  (4.25) 
for  a  =  0,  are  multivalued  functions.  Their  principal  branches,  as 
demonstrat  i  in  Fig.  I  and  accompanying  discussion,  are  given  by  (2.5)  and 
(2.4),  respectively,  which  accounts  for  the  adoption  of  the  convention  (2.6), 
and  which  for  our  purpose  we  rewrite  here  as  follows: 

(-w)c  =  ^Cg-i^Cc^  e  =  sgn  { Iinw  1  ,  (5.7) 


where  c  Is  an  arbitrary  exponent.  We  propose  to  apply  this  rule  to  compute 


the  powers  (-i£)  a  +  and  (i£)  3~  which  appear  in  (5.1)  and  (5.2).  Thus,  when 

Imfis)  >  0  in  (5.1),  we  have  e  =  1  and  (~i£)  3  +  ~  (-i£)  3+  according  to  (5.6), 

whence 

(-U)-a+  =  (i?)_a+ei1Ta+;  el7Ta+  =  el7T(a/2  +  i?o>  =  ^  e~^o  ,  (5.8) 

and  when  Lm{i£}  <  0  in  (5.1),  we  have  simply  (-!£)  a+  ~  (i£)~a+.  Similarly, 
when  Ira{iE;}  >  0  in  (5.2),  we  have  (i£)  a~  ~  (i£)-a“,  and  when  Im{iO  <  0, 

we  have  e  =  -1  and  (i5)*"3~  ~  (-i£)“a“;  i.e., 

(-i?)"a~  =  (i?)"3-  e~i1,a“;  e-i7Ta“  =  e-i"(a/2  -  i?0)  =  i~a  e-7TB0  #  (5.9) 

In  (5.8)  and  (5.9),  we  still  have  to  compute  the  powers  (i£)  a±,  where 

a+  =  t  ±i3  .  Thus,  we  have 
l  o 

(i?)-3-*-  =  e~a+An(i5)  _  fi“( ot/2  +  iB0)in(iO  _  ^^-a/2  g ^ 0O / 2  -  i30ZnC.  ( 5.  iq) 

(  i  r)  ^ —  =  g~a— ^n(i^)  —  (  ®/ 2  —  i 80) £n(i £)  _  ct/2  tt80/2  +  iSQlJ.n'^  (5.11) 

The  logarithmic  terms  in  the  last  exponentials  in  (5.10)  and  (5.11),  when 
combined  eventually  with  e-i^/2,  are  asymptotically  ( £  +  »)  negligible  and 
they  are  frequently  omitted.  Here,  however,  we  retain  the  logarithmic  terms 
for  completness  sake. 

Hence,  combining  the  results  of  the  two  preceding  paragraphs  with  our 
basic  leading  term  asymptotic  expansions  (5.1)  and  (5.2),  we  may  now  tabulate 
the  asymptotic  values  of  Z+(£)  and  Z_(  £)  for  Im{i£}  >  0  and  Im{i£}  <  0,  which 
is  equivalent  to  saying  z  >  0  and  z  <  0,  respectively.  Thus,  we  obtain  the 


four  expressions: 


ImUC)  >  0;  Z+(C)  ~ 


i(1  a)  r(a  +  )[iae  n3°/2J(iU  r-/2"i  3<>*n: 


Im(i C }  <  0;  Z+(C)  -  i(1  a)  T(a  + )  [e  "  2  )  ( i  l)  ^V’'2 


IrntiC}  >  0;  Z  (C)  -  -i(1  a)  Z(a  )  [e"*  3°'' 2  ]  ( i  l) ' V  V  1 2+1  ">ln 


ImUU  <  0;  Z 


.(C)  ~-i(1"a)  r(a  )[i  V3,B^;i](i()  o.It) 


in  which  the  first  two  represent  progressive  waves  and  the  last  two  regressive 
waves.  The  given  limiting  forms  are  valid  for  all  positive  values  of  the 
parameter  a,  0  <  a  <  ®.  For  Budden's  problem  corresponding  to  (5.3)  and 
(5. A),  we  put  a  =  0  everywhere  except  in  the  normalization  factors  were  we  put 
a  =  2  in  accordance  with  (A. 23). 


To  conclude  this  Section,  we  extract  from  the  preceeding  equations  the 
spatial  dependence  of  the  wave  functions  of  unit  amplitude,  namely: 


(C)  ~  (iC)'a/2eU/2  i3o*"-'  =  (ir)  a+eir’/2  ;  (5.1b) 


(?)  ~  (iC)~a/2e'i  ’/2+i3otn’  =  (iC)-a-e-i  CJ1  ;  (5.17) 


which  equally  well  may  be  expressed  as  functions  of  the  real  independent 
variable,  -«  <  z  <  »,  albeit  no  longer  of  un i t  amplitude;  that  is: 


*+(z)  ~(2iB2)  3+  |  z  |  a+ei32Z  32 


♦  (z)  ~( 21 0O)  a-  |zj  a~e 


(5.1b) 


c  s . is> 


-V  V/.*.',').1  V  V  V  Vs." W  V  ‘ 


which  are  Che  functions  that  supplant  the  forms  given  in  Eqs.  (7)  in  Ref.  1. 

It  is  important  to  reiterate,  for  both  pairs  of  equations,  that  a+  =  a/2  t  i30 
5  =  202(z  +  i P) ,  and  %  =  2 B 2 1 z I »  where  arg  5  =  arg  02  2 


Finally,  to  conform  with  the  nomenclature  introduced  in  Ref.  1,  we 
present  in  Fig.  4  a  diagram  indicating  the  directions  of  propagation  of  the 
four  asymptotic  leading  terms  for  Z+(5)  and  Im{i£}  ^  0,  which  are  associated 
with  the  uni t  amplitude  wave  functions  (5.16)  and  (5.17).  That  is,  in  the 
spirit  of  Eqs.  (7)  of  Ref.  1,  we  have 


iL(C)  -  (iO 


-(a/2-i0o)  -U/2  . 

"  f 


(5.20a) 


»,(5) 


(iO 


-( a/2+i 0O)  15/2  . 

0  e  ; 


(5.20b) 


*3(5)  ~  (U)  <a/2+l6o)e^/2 


(5.20c) 


MO  -  ( i  5 ) 


i?)-(a/2-i0o)e-i5/2 


( 5. 20d ) 


It  is  seen  from  Eqs.  (5.20)  that,  when  a  >  0,  there  is  introduced  a  decaying 
envelope  of  the  form  ?~ai/2  =  (28  |z|)  2 ,  or  conversely,  it  can  be  inter¬ 

preted  (asymptotically)  as  an  enhanced  swelling  of  the  wave  form  toward  the 

-2 

resonance,  which  we  associate  with  term  5  in  the  effective  potential  (3.4). 
It  is  worth  noting  that  the  swelling  disappears  when  a  =  0,  which  is  Budden's 
problem.  In  addition,  it  is  pertinent  at  this  juncture  to  emphasize  that  the 
exponentials  in  (5.20),  exp{±U/2}  =  exp{±i02z  +  02F},  contain  the  constant 
attenuation  factors  expf+B-jT}  =  exp{+h/2}  in  accordance  with  the  definitions 


The  connection  that  exists  between  the  unit  amplitude  wave  functions 
listed  in  (3.20)  and  the  actual  asymptotic  solutions  Z+(r)  that  we  derived  lor 
lm{iC)  ^  0  is  as  follows:  the  progressive  waves  and  '3 C  ‘w )  are  contained 

in  the  functions  Z+( £)  defined  by  (3.12)  and  (5.13),  respectively;  similarly, 
the  regressive  waves  $^(£)  and  ^(d)  are  contained  in  the  functions  Z_(r) 
given  by  (5.14)  and  (5.15),  respectively.  The  multiplying  coefficients  in 
these  equations,  as  we  shall  see,  are  of  the  essence  in  the  computation  of  the 
transmission  and  reflection  coefficients  that  we  undertake  at  once. 
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b.  Computation  of  Transmission  and  Reflection  Coefficients 

The  computation  of  the  reflection  and  transmission  coefficients  associated 
with  the  leading  term  asymptotic  forms  for  the  solutions  Z  +  (£)  of  our 
generalized  equation  (3.1)  is  conveniently  carried  out  by  resorting  to  the 
integral  representations  that  we  derived  in  Section  4,  and  by  applying  the 
method  of  analytic  continuation.  This  means,  with  reference  to  Fig.  4,  that 
we  first  set  up  asymptotic  transmitted  waves  of  unit  amplitude,  ^2( ?)  lor 
z  >  0  or  4>4(C)  for  z  <  0  (cf.  remarks  preceeding  (5.12)),  and  then  we  proceed 
to  seek  their  analytic  continuation  as  z  changes  sign.  In  what  follows  we 
shall  assume  temporarily  that  £  =  202Z  is  rigorously  real  and  covers  the 
infinite  domain,  <  £  <  °°.  Accordingly,  the  four  half-ray  rectilinear  paths 
of  integration  defined  by  (4.5)  and  (4.13)  may  be  drawn  for  convenience  as 
strictly  vertical  in  the  complex  t  plane,  as  shown  graphically  in  Fig.  3. 

First  we  consider  the  left  to  right  case  of  Fig.  4  and  then  making  use  of 
(4.5),  we  draw  the  diagram  of  Fig.  5,  which  illustrates  the  two  rectilinear 
paths  of  integration  corresponding  to  Z+(  E,  >  0)  with  $  =  0  and  Z+(  £  <  0)  with 
<{>  =  it,  and  indicating  that  their  unit  amplitude  wave  functions  are  and 

<t>  3  ( C ) ,  respectively.  Thus,  as  stated  in  the  preceding  paragraph,  we  identify 
the  linear  path  of  integration  for  Z+( £  >  0)  as  the  transmitted  wave  for 
z  >  U.  The  analytic  continuation  of  this  solution  for  z  <  0,  which  is  the 
incident  wave  Z+( £  <  0),  is  obtained  by  rotating  the  half-rav  path  of 
integration  clockwise  about  the  branch  point  at  t  =  1/2,  as  indicated  by  the 
(directed)  semicircle  on  Fig.  5.  This  procedure  is  dictated  bv  the  unit 


circle  in 

Fig.  3(a),  which  shows  that 

t  lie 

t  rans i t i on 

from  t> 

-  arg  \  =  ■)  to 

<t>  =  arg  z 

=  tt  is  achieved  by  means  of 

the 

thick  semi c 

i re le  t 

o  tile  right  ot 

vertical. 

Further,  this  procedure  imp 

1  i  e  s 

an  exon  rs  i 

on  into 

another  sheet 
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of  Che  Riemann  surface,  which  we  indicate  in  both  Figs.  3(a)  and  5  by 
drawing  the  path  for  Z+( 5  >  0)  as  full  line,  and  the  path  for  Z+( £  <  0)  as 
dashed  line.  The  transition  from  the  incident  wave  3 (  ^ )  to  the  transmitted 
wave  4>2(?)  is  effected  by  means  of  an  infinitesimal  semicircle  traversed  in 
the  counterclockwise  direction  about  the  branch  point  at  t  =  1/2,  which 
introduces  the  factor  exp(Ur(a/2  +  i 3a ) }  that  we  identify  below  as  the 
reflection  coefficient. 


Hence,  Z+( £)  is  its  own  analytic  continuation  as  we  change  the  sign  of  z, 
and  in  the  left  to  right  case  of  Fig.  A,  there  is  no  reflection,  a  result 
also  found  in  the  Budden  problem  (a  =  0).  Thus,  the  reflection  coefficient 
for  <j>3(£)  going  into  >4(0  vanishes  identically  for  all  a,  =  0.  To 

compute  the  transmission  coefficient  T32  from  d) 3 (  \)  into  $4(3)  quite 
generally;  that  is,  in  the  presence  of  dissipation,  we  first  recall  from  the 
discussion  that  follows  (5.20)  that  the  unit  amplitude  wave  functions  (3. In) 
and  (5.17)  now  exhibit  the  attenuation  factors  e-r'/2j  and  then  we  take  the 
ratio  of  coefficient  of  $+(£)  =  $2(11)  in  (5.12),  the  transmitted  wave,  to  the 
coefficient  of  $+(  £)  =  $3(0  in  (5.13).  Hence,  taking  due  account  of  the 
attenuation  factors  just  mentioned,  we  obtain 


(  l-a)  n,  \ r  .  ot  -it  3  /  2  1  -n/2  1  -n.3  /  2 

i  1 (a  )  [ 1  e  o  e  i  e  o 


(l-a)  r  tt 3  / 2  1  -n/2 

T(a+)Le  o  ,  e 


"3  /'2 

l*  (  ) 


rVnBo  • 


( t>.  1 ) 


which  is  completely  equivalent  to  the  factor  already  given,  and  which  shows 

1  i  TT  ■  1/  * 

explicitly  the  periodic  dependence  on  t  through  the  phase  factor  i  -  e 
Thus,  for  Budden's  problem  (1  =  U),  we  recover  the  well-known  result 


(  h .  2  ) 


and  for  the  self-adjoint  analogue  or  budden's  equation  (a  =  1),  we  obtain 

T32  =  ie_1T .  (6.3) 

Notice  that  all  three  values  of  the  reflection  coefficient  T32  given  above  are 
contained  in  the  factor  exp{iw(a/2  +  i ) }  that  arose  from  the  infinitesimal 
semicircle  about  t  =  1/2  when  traversed  in  the  counterclockwise  direction. 

Next,  we  consider  the  right  to  left  case  of  Fig.  4,  i.e.,  the  incident 
wave  <i> x ( 5 )  approaches  resonance  from  the  tunneling  side  (2  >  0)  and,  as  we 
shall  see,  gives  rise  to  a  transmitted  wave  414(C)  for  z  <  0  and  a  reflected 
wave  -(1  -  e2  ia+)  42(C)  for  2  >  0.  Thus,  proceeding  as  before,  we  make  use 
of  (4.5)  and  (4.13)  to  draw  the  diagra  of  Fig.  6  illustrating  the  three 
rectilinear  paths  of  integration  corresponding  asymptotically  to  the  Incident 
wave  Z_(C  >  0)  for  4  =  0,  the  transmitted  wave  Z_(  C  <  0)  for  4>  =  tt ,  and  the 
ref  lected  wave  -( 1  -  e^1Tia+)  Z+(  r  >  0)  for  4  =  0  and  4  =  2 tt.  We  have 
identified  the  rectilinear  path  of  integraton  for  Z_( C  <  0)  as  representing 
the  transmitted  wave  for  z  <  0.  The  analytic  continuation  of  this  solution 
for  z  >  0  is  obtained  by  rotating  the  half-ray  path  of  integration 
counterc lockwlse  about  the  branch  point  at  t  =  -1/2  as  indicated  by  the 
(directed)  semicircle  in  Fig.  5.  This  procedure  is  dictated  by  the  unit 
circle  in  Fig.  3(b)  which  shows  that  the  transition  from  4  =  arg  C  =  "  to 
4  =  arg  C  -  0  is  achieved  by  means  of  the  thick  semicircle  to  the  right  of  the 
vertical.  In  so  doing,  however,  we  cannot  ignore  the  presence  of  the  branch 
point  at  t  =  L/2,  with  the  result  that  the  proposed  deformation  of  the  path 
of  integration  yields  the  incident  wave  Z_(  c,  >  U)  from  (a.l  i)  with  t>  =  U,  plus 
the  reflected  wave  from  (4.5)  with  4>  =  U  and  4  =  2tt,  which  we  write  as 
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’  t  i  a  , 

-(  1  -  e~  )  Z+(s  >  u).  Here  a+  =  a/ 2  +  i  80  is  the  exponent  associated  with 


the  branch  point  at  t  =  1/2.  We  observe  that  the  contour  of  integration  for 
the  reflected  wave  is  represented  in  Fig.  6  by  a  solid  halt-line  transversed 
downwards,  an  infinitesimal  circle  around  the  branch  point  at  t  =  1/2 
traversed  in  the  coun t e r c loc kwi se  direction,  and  a  dashed  half-line  traversed 
upwards  in  another  sheet  of  the  Riemann  surface.  We  also  observe  that  the 
transition  from  the  incident  wave  function  Z_( r  >  0)  to  the  transmitted  wave 
function  Z_(  r,  <  0)  is  effected  by  means  of  an  infinitesimal  semicircle 
traversed  in  the  clockwise  direction  about  the  branch  point  at  t  =  -1/2,  and 
giving  rise  to  the  factor  exp { -i ^ ( a/ 2  -  iS0)}.  Based  on  the  results  of  the 
preceeding  paragraph,  we  identify  this  factor  with  Tj/,,  the  transmission 
coefficient  from  t>  |  (  r, )  into  i>4(%);  that  is, 


.-"3u 


(6.4) 


which  is  the  complex  conjugate  ot  (6.1)  m  the  absence  of  dissipation,  when 
2,,  i real.  As  bef  >re,  we  consider  two  cases.  For  Budden's  problem  we  put 
i  =  i'i  (  h .  )  to  obt  ain 


(6.3) 


snub  is  Budden's  result  and  is  identical  to  (6.2).  For  the  se  It -adjoint 
anal  igue  >f  Budden's  equation,  we  pur  i  =  1  l:i  <6. a)  to  yield 


(  n.  b) 


1 1  i  s  the  rnrii  !"x  i  in  hi,  ate  o  t  (6.3)  w :  i  e  a  2  is  real. 


We  have  shown  chat  the  computation  of  the  transmission  coefficients  T32 


and  T14  emerges  quite  simply  from  the  contours  of  integration  and  does  not 
necessitate  a  detailed  knowledge  of  the  asymptotic  solutions  given  by  the 
integral  representations.  The  same  statement,  however,  does  not  apply  to  the 
computation  of  the  reflection  coefficient  Rj2«  In  fact,  we  must  proceed  as  we 
did  originally  with  the  computation  of  the  transmission  coefficient  (b.l)  in 
which  we  included  dissipation  and,  making  use  of  the  conclusions  following 
(5.20),  we  took  cognizance  of  the  fact  that  the  unit  amplitude  wave  functions 
(5.16)  and  (5.17)  now  exhibit  the  constant  attenuation  factors  eFri/2^ 
respectively.  When  computing  transmission  coefficients  by  taking  ratios  of 
the  pertinent  coefficients  in  (5. 13)  -  (5.15),  we  find  that  these  factors  play 
no  role  because  they  cancel  each  other.  However,  this  is  no  longer  the  case 
when  computing  the  reflection  coefficient  Rj2*  Thus,  we  must  include  in 
addition  to  the  factor  -(1  -  e^iria+)  shown  in  Fig.  6,  the  ratio  of  the 
coefficient  of  i>+(^)  =  i>2C  ? )  in  (5.12),  to  the  coefficient  of  d>_(  5 )  =  4>i(0 
in  (5.14),  and  we  must  take  due  account  of  the  constant  factors  e?r|//^,  which 
arise  in  the  presence  of  dissipation.  In  this  manner,  we  obtain 


,  .  i(l_a)r(a  )  [iae_,r8o/2]e-n/2 
R12  =  -(1  -  e2"13*)  - - - - - - 

-i(L_a)r(a_)  [e“nBo/2]e-n/2 


T  ( a/  2  +i  8  ) 

- 1"1  — : 

T( 3/ 2-i 8  ) 
o 

which  is  formally  identical  to  Eq .  (10)  of  Ref.  1,  except  for  the  fact  that 
here  n  =  232r,  whereas  previously  we  wrote  v  instead  of  T.  From  (6.7),  we  can 
deduce  the  special  case  corresponding  to  the  self-adjoint  analogue  of  Budden's 


i  -  e 


ilT3  —  2lT  R  ,  -n 


!°  ]« 


(0  <  a  <  00 ) 


(6.7) 


equation  by  putting  a  =  1;  that  is, 


-  3  i- 


tn  which  ]  R  2  i  >  1  unless  n  is  extremely  large.  Finally,  tor  Budden's  case 
(a  =  (.1),  we  observe  that  we  rnav  put  a  =  2,  or  else  a  =  U,  in  all  terms  ot 
( o • 7 )  to  obtain  the  equivalent  torrns^ 


-2n3o 

in  accordance  with  Budden's  original  result,  but  now  showing  the  explicit  role 
of  dissipation  through  the  complex  parameters  Bq  and  n.  Notice,  furthermore, 
that  in  (b.9)  we  have  | R  j  2 1  x  1  under  all  circumstances. 


e  n  ,  (b.9) 


2(1  +  18  ) 

R  = - -°-  [1  - 

r(l-iSo) 


*P(  i  3  ) 

]e“n  4-  [1  "  e 

r(-iBo) 


The  most  unexpected  result  of  the  present  investigation  is  the  discovery, 
embodied  in  r.q .  (b.7),  that  the  reflection  coefficient  R;;  exhibits  an 
anomalous  behavior  as  a  function  ot  a.  To  study  this  effect  more  thoroughly 
and  to  simplify  the  computations,  we  confine  our  attention  in  further  remarks 
to  the  case  of  no  dissipation;  this  means  that  in  (b.7)  8^  =  3  =  8^(23^) 

is  real  and  1  =  23.,  T  =  o.  Hence,  we  deduce  that 


I  K  j  ?  j  ^  =  1  +  e  -  2cos(aw)e  “1T^o  , 


(b.  10) 


which  shows  that  the  square  of  the  modulus  of  the  reflection  coefficient, 

K 1  2  ‘  t  oscillates  periodically  as  a  is  varied,  U  <_  a  <  00 .  The  smallest 

.  -2^8 

reflection  occurs  for  1  -  2N,  N  =  0,1,2,...,  which  is  min  [Kjuj  =  ]  -  o , 

whereas  the  maximum  reflection  takes  place  when  a  =  2N  +  l ,  N  =  0,1,2,...,  and 

-2*3 

att  iins  the  value  max|K[yl  =  1  +  e  o,  which  is  greater  than  unity.  This 
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enhancement  in  the  reflection  coefficient  indicates  that  there  exist  bands  of 
values  of  a  within  which  the  resonance  spontaneously  emits  waves.  The  onset 
for  stimulated  emission  is  obtained  from  the  condition 

!  T 1 4  I  2  +  I  Ri2  I  2  =  0  ,  (6.U) 

from  which,  making  use  of  (b.4)  and  (b.lU),  we  deduce  the  equation 

2cos( am)  =  1  +  e”2wB°  ,  (b. 12) 

whose  roots  determine  the  values  of  a  at  the  edges  of  the  emission  bands.  We 
note  further  that 

I  R 1  2  (  I  2  =  IR12^'2  +  4  sln2(  aw/ 2)e  ^1TS°  >  (  b.  13) 

which  indicates  that  except  for  the  values  a  =  2N,  N  =  0,1,2,..., 
corresponding  to  the  minimum  reflection,  the  result  of  putting  a  *  2N  is  to 
enhance  the  reflection  coefficient  above  the  Budden  value,  |Rj2(u)i2* 

In  summary,  the  presence  of  the  first  derivative  term  a<i>'(z),  0  <_  a  <  <*>, 
in  the  generalized  resonance-tunneling  equation  (2.1)  that  we  have  constructed 
from  the  well-known  Budden's  equation,  can  yield  both  absorption  as  well  as 
stimulated  emission  of  waves,  depending  on  the  periodic  band  conditions  on  a 
that  we  established  here.  It  would  be  interesting  and  challenging  to  find 
examples  of  physical  systems  where  such  a  behavior  occurs. 

This  work  has  been  supported  by  the  Office  of  Naval  Research. 
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